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ABSTRACT 

Aims. The origin of the observed variability of the gas-phase D/H ratio in the local interstellar medium is still debated, and in 
particular the role of deuterium depletion onto dust grains. Here we extend the study of the relationship between deuterium and 
titanium, a refractory species and tracer of elemental depletion, and explore other relationships. 

Methods. We have acquired high resolution spectra for nine early-type stars using the VLT/UVES spectrograph, and detected the 
absorption lines of interstellar Till. Using a weighted orthogonal distance regression (ODR) code and a special method to treat non 
symmetric errors, we compare the Till columns with the corresponding HI, DI and also 01 columns. In parallel we perform the same 
comparisons for available Fell data. 

Results. We find a significant correlation between Till/Hl and D/H in our data set, and, when combined with published results, we 
confirm and better constrain the previously established trends and extend the trends to low HI columns. 

We exclude uncertainties in HI and OI columns as the main contributor to the derived metals-deuterium correlations by showing that 
the Till/HI ratio is positively correlated with DI/OI. We find a similar correlation between Fell/HI and DI/OI. 

The Till gradients are similar or slightly smaller than for Fell, while one would expect larger variations on the basis of the higher 
condensation temperature of titanium. However we argue that ionisation effects introduce biases that affect iron and not titanium and 
may explain the gradient similarity. 

We find a less significant negative correlation between the Till/DI ratio and the hydrogen column, possibly a sign of different evapo- 
ration of D and metals according to the cloud properties. More Till absorption data along very low H column lines-of-sight would be 
useful to improve the correlation statistics. 

Key words. Interstellar Medium - Deuterium Abundance - 
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j> ' 1 . Introduction (LOS) with log N(HI) < 20.5 cm"2, D/H is significantly lower, 

• --j ^7 ppm. 

y\. The measurement of the present-day ratio of deuterium (D) to . j- ^ • ^-irvrv/ix j t - i . i /'-irvrv^\ / 

, , • , ^ , , • . • According to Drame (2004) and Lmsky et al. (2006) (see 

H hydrogen (H) in the Galaxy places important constraints on Its , , • i '<hnn\ j • u r%/ii 1- ■ 

Cn , .'^ , ' J , , , , r ■ , r also Linsky 2007), the wide range in the gas-phase D/H ratios is 

. . chemical evolution and has been the goal of intense work, from , ^ ' , ^ ■ j w . • . . n 

, _ . ,„ ■ , ,,f,r^ / due to different amounts of deuterium depletion onto interstellar 

the Copernicus era (Rogerson and York, 1973) to HST (eg. . ™ . , .u- • . . • c 

T- 1 Ti. « « r.r, /T 1 • 1 -.ni^r^ o 1 1 grains. Theoretical support for this interpretation IS the preferen- 

Linsky, 1998), IMAPS (Jenkins et al., 1999, Sonneborn et al., ^, . f. . • . • . • • , 

jr-irnr-/ »j i »j • j tial adsorption of deutcnum onto interstcllar grains, in particulaT 

2000) and FUSE (e.g. Moos et al., 2002). Measurements inside ,, i„oo • onn/i\ /-.u i . 

, , ,. PAHs (Jura, 1982, Draine, 2004). Observational support comes 

the Local Bubble (LB), the volume of tenuous gas surrounding „ . . , ^ j o- j i 

, „ , „ , , ,„„o , ■ , ,^^r^ T 11 from the coiTelations between D/H and both Fe and Si deple- 

the Sun (e.g. Snowdenetal., 1998, Welsh etal., 1999, Lallement . „ , . , . / 

1 1 • 1 -1 1 .c ,1 1 tions, as wcll as the coiTelations with the rotational temperatures 

et al., 2003) appear to be consistent with a single value for D/H nr r i- r • • i i /t • i . i 

1 It -1 J D »i ^1- • 1 of H2 for lines-of-sight posscssing molecular gas (Linsky et al., 

(Moosetal., 2002, Hebrard& Moos, 2003, Oliveira etal., 2003), ^^.^.A t.. , a .a ^ u \ ■ a u .u 

, , , , 1 . 1 , o , 2006). The D/H trends quoted above may be explained by the 

but other measurements (e.g. Jenkins et al., 1999, Sonneborn et . i ' • u • r u • /u r.u 

,1 -^nnn h„..,»c »t ,i -^nn^ »t ,i ^hha h^k.-,.h depletion hypothesis if one hypothesizes a recent heating of the 



al., 2000, Hoopes et al., 2003, Friedman et al., 2006, Hebrard 
et al., 2005, Oliveira & Hebrard, 2006) suggest variations of the 
interstellar D/H ratio at larger distances, i.e. beyond ^ 150-200 



LB and especially recent shocks in star-forming regions at its 

periphery. This interpretation has important consequences and 

„ -I •,• J 1 •, .cirn implies that the total D/H ratio including deuterium in the gas 

parsecs. Quantitatively, within and at the periphery of the LB a- .i. a .^ . . n- j u • u 

/I xt/tt/s 1^^-, 9n 1 r^/TT • • 1 1 and in thc dust IS at Icast 23 parts pcr million of hydrogcn, wluch 

(log N(HI) < 19.3 cm"2), the D/H ratio is of the order of 15 . „ . j , ^ i u • i 1 

, L. , ^,.1 1 ^o„^s 111. J- • 1 poses a challenge for models of galactic chemical evolution, 
ppm (see Figure 1 of Linsky et al., 2006), while lines-of -sight 

measurements with higher N(HI) show a strong variability, with On the other hand, Oliveira & Hebrard (2006) show that high 

maximum values of ^ 23 ppm. For a number of lines-of-sight D/H values are also found along distant lines-of-sight with asso- 
ciated high interstellar columns. At high N(HI) the selection bias 

against high D/H values (due to dominance of dense, cold, and 

Send ojfprint requests to: R. Lallement thus Strongly depleted regions) should not lead to any variability 
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(Oliveira et al., 2006). Moreover, D/O measurements do not fully 
support this depletion scenario, since the DI/OI ratio is found to 
be more constant than DI/HI (Hebrard et al., 2005). Interestingly 
the LB D/H ratio derived from the D and O measurements is 
only of the order of 13 ppm (Hebrard & Moos 2003). These au- 
thors suggest that the uncertainties in the measurement of HI in 
the data play an important role in the subsequently derived re- 
lationships. Two quotients with the same denominator remain 
naturally correlated if there are measurement errors on this de- 
nominator. Fell/HI vs. DI/HI or Fell/OI vs. DI/OI relationships 
could result from this effect. 

Infall of external unprocessed (or less processed) gas also 
creates D/H variability, and indeed a continuous infall is a gen- 
eral property in the Galaxy (Romano et al., 2006). Support for 
the existence of a significant infall into the Sun's neighborhood 
comes from D, O, N, and ^He abundance values found in the lo- 
cal interstellar cloud that surrounds the Sun. Geiss et al. (2002) 
interpret these abundances as being due to the local mixture of 
processed galactic gas and very weakly processed (Magellanic 
type) gas of intergalactic origin having fallen onto the disk. Such 
an infall could have occurred recently, leading to a present-day 
incomplete mixing. On the other hand, the wide D/H range im- 
plies a non realistic infall strength and a very weak mixing after 
the infall, as argued by Linsky et al. (2006). 

It is therefore mandatory to investigate in more detail the 
depletion mechanisms and the associated deuterium behaviour. 
Prochaska et al. (2005) made an important step by showing a 
positive correlation between the abundance of titanium, a re- 
fractory species, and that of D/H. Here we present results we 
have obtained in the southern hemisphere with the Ultraviolet 
and Visual Echelle Spectrograph (U VES) at the ESO Very Large 
Telescope (VLT), and we add nine new Till measurements to- 
wards D/H targets. In parallel to our observations and analysis, 
Ellison et al. (2007) also gathered and analysed new VLT-UVES 
and Keck Observatory HIRES data. We have combined our re- 
sults with their results and those from Prochaska et al. For tar- 
gets in common, we have compared the two results. We have 
extended the correlative study, in particular to lower column- 
densities. We have used sophisticated methods to treat asymmet- 
ric errors in both correlated quantities. 

Titanium is a particularly well suited element for the study 
of a potential correlation between depletion and deuterium abun- 
dance. It has a high depletion which is Unked to its high con- 
densation temperature. Till is optically accessible in the near 
ultraviolet via the AA 3229.199,3242.994,3383.768 A absorp- 
tion lines (wavelengths from Morton, 1991), and interstellar ti- 
tanium has been observed in absorption since the 1970's (e.g. 
Wallerstein and Goldsmith, 1974, Stokes, 1978, Hobbs, 1984). 
Moreover, as already argued by these authors, Ti II is the dom- 
inant ionization state in HI regions due to the nearly exact co- 
incidence of the ionization potentials of Till and HI. Also, 
Ti II is coupled to HI via charge-exchange reactions (Field and 
Steigman, 1971). The observed close N(Till) vs. N(HI) rela- 
tionship, in contrast with the more dispersed N(NaI)/N(HI) or 
N(CaII)/N(HI) relations (see e.g. Figure 7 of Welsh et al., 1997), 
results from those physical properties. 

The detection of the interstellar Till Unes requires high spec- 
tral resolution (R > 50,000) and appropriate optical systems 
and detectors adapted to the near ultraviolet wavelength range. 
Prochaska et al. (2005) used HIRES at the Keck Observatory 
and observed absorption towards 7 northern hemisphere stars. 
Due to this limited number of targets, it is important to ex- 
tend their study to other D/H targets. Moreover, the observed 
trend is strongly dependent on the measurement towards Feige 



110, and indeed Hebrard et al. (2005) have shown that the col- 
umn of HI is probably significantly underestimated for this star. 
They also quote the contradicting case of y^Vel, a line-of-sight 
with high D/H and low titanium abundance. Finally, while they 
covered a range of logN(HI) from 19.85 and 21.05 cm"^, it is 
important to extend the study to lines-of-sight with lower HI 
colutims, because there is a strong D/H variabihty for logN(H) 
< 20 cm"^ and no titanium column density data are currently 
available within this range. 

Section 2 describes the observations, data analysis and Till 
column determination. In section 3 we perform Unear compari- 
son fits between Till and other species, using our data in com- 
bination with other pubhshed data. We perform corresponding 
correlative studies for available Fell data. Section 4 discusses 
the results. 



2. Observations and absorption line analysis 

The spectroscopic data were obtained with the ESO VLT-UT2 
telescope and the UVES spectrograph (Dekker et al., 2000) un- 
der the observing program 077.C-0547A. We used the standard 
dichroic 346-5 80nm setting (DICl) and the image sheer 1 which 
corresponds to R ^ 60,000 at 3400 A. This setup included the AA 
3229,3242,3384 A Ti II lines with the blue arm. The spectra were 
recorded in service mode. 

The targets were seven bright early-type stars and two hot 
subdwarfs. They were selected from the list of D/H targets stars 
recently compiled by Linsky et al., (2006). Multiple short expo- 
sures were obtained for each of the bright stars and two long ex- 
posures for the subdwarfs. Details of the observations are listed 
in Table 1 . The individual stellar spectra, processed through the 
UVES pipeline, have been co-added. 

The spectral analysis and the correlative studies were per- 
formed in the frame of the IGOR/WaveMetrics software, using 
both built-in and specially developed routines. In the present 
work we are interested in the total Ti II column-density values (a 
cloud-by-cloud analysis is beyond its scope), and since we are in 
the optically thin linear regime measurement of the total equiva- 
lent width is sufficient. For this purpose, stellar continua around 
the strongest transitions AA 3383.768, 3241.994 A were fit with 
a multi-order polynomial. The noise level was automatically cal- 
culated during the stellar continuum fit and assumed to prevail at 
the absorption Une locations. The resulting signal to noise ratios 
around the 3 3 84 A transition are given in Table 1 . The normalised 
spectra corresponding to the two strongest transitions are dis- 
played in Figure 1. Note that these targets stars correspond to a 
two orders of magnitude variation of the HI column-density. 

After division by the fitted continuum, equivalent widths are 
simply deduced from the normalised spectra. Errors on the total 
Ti II equivalent width are estimated on the basis of the noise level 
and the number of pixels defining the absorption. 

Two of the target spectra (6 Car, Pup) are afifected by some 
instrumental interference fringes. For these stars the stellar con- 
tinua were fit with the sum of a multi-order polynomial and a 
sinusoidal function with an adjustable period, phase and amph- 
tude. The correction is illustrated in Figure 2 for 6 Car which 
corresponds to the maximum fringe amplitude. The diff'erence in 
derived equivalent widths (or column densities) with and with- 
out the fringe correction is generally less than the error resulting 
from the noise. 

The column-densities of Till are deduced here from the 

strongest 3384 A transition. On the basis of experimental data 
Pickering et al (2001) have revised the oscillator strength of this 
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Fig. 1. Normalised spectra around the Till 3384 A (solid lines) and 3242 A (dashed lines) transitions for our target stars 

. The abscissa is the heliocentric velocity in km/s. 



transition at f=0.358, some 5 percent higher than the Morton 
(1991) value f=0.3401. In order to be consistent with the previ- 
ous analyses of Prochaska et al. and Ellison et al. we have used 
this revised value. Note that the results of this work and previ- 
ous ones on the existence and quality of correlations remain un- 
changed if this oscillator strength is revised again. Only resulting 
gradients would have to be scaled according to the changes. 

For the AA 3242 A transition we first used the Morton (1991) 
atomic data, then the revised values quoted by Prochaska et al 
(2005). Their atomic data obtained from Morton (2003, private 
communication) are significantly different from the older ones, 
with the oscillator strength increased from f=0.183 to 0.232. We 
could check from the simultaneous fits of the two lines that those 
revised values are in better agreement with data than the Morton 
(1991) values and we show only the second results. 

We have also performed line-fitting, using the program de- 
scribed in Sfeir et al. (1999), i.e., the residual intensity absorp- 
tion profiles were fit to combinations of Voigt profiles, convolved 
by the UVES instrumental function. We have limited the temper- 
ature to the 100-15, OOOK range, which corresponds to realistic 
values in the local ISM, and we have simply chosen the mini- 
mum number of clouds allowing a good visual adjustment to the 
data. In parallel we performed a profile fitting analysis with the 
Owens. f software of Lemoine et al. (2002) for one of the targets 
and the results are consistent. The cloud-by-cloud structure in 
those line-fitting results should not be considered as the best so- 
lutions, since they should also be compared in future work with 
ithe structure derived from other elements. However, they pro- 
vide a consistency check. Simultaneous fitting of the two transi- 
tions are shown in Figures 3 and 4. We have not included aVir, 



jSCen and aCru, for which the column is small and the detection 
of the weakest Une is of insufficient quality and does not add to 
the 3384A detection. 

Although we do not specifically use the column densities or 
the cloud-by-cloud structure resulting from this exercise, fitting 
the two 3384 and 3242 A lines, as said above, has provided a 
consistency check: Lines other than Till and (or) non-random 
noise features are revealed by this procedure. 

While performing this exercise we noticed the following 
departures from the model absorptions. For LS1274 the AA 
3241.994 A region is affected by an additional absorption that 
has no counterpart at 3384A not at 3229 A (the absorption is 
strong enough for this star to be detectable at this wavelength). 
Looking at the entire spectrum, we found that at variance with 
the other targets and especially the second subdwarf JL9, the 
spectrum is characterized by a very large number of narrow fea- 
tures. We could not determine whether these narrow lines have a 
circumsteUar or stellar origin (LS 1274 is cooler than JL9) origin, 
and their study is beyond the scope of this paper. For this star we 
have calculated the column-density without including this addi- 
tional absorption that is not seen in the stronger transition and 
is indicated in Figure 1. A second, although smaller, anomaly is 
found for aCru, for which there is an absorption feature seen at 
3342 A and not for the strongest transition. The noise level how- 
ever is such that a non-random noise feature may be responsible. 
After examination of the two transitions and of these line fit re- 
sults we have accordingly increased the error bars on the equiva- 
lent widths and resulting columns to take into account such non 
random noise features. 
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Fig. 6. Titanium vs. deuterium colunms. The star /zCol is clearly 
an "outlier". 



Fig. 2. Illustration of fringe correction to the data: the extreme 
case of dCai. 

We also have searched for a possible contamination by a stel- 
lar line by shifting each spectrum to the stellar reference frame 
and overplotting the resulting shifted spectra. No particular ab- 
sorption at a given wavelength is revealed by such a method. 
The two subdwarfs however are excluded because stellar radial 
velocities for those targets are uncertain. 

For all targets except fiCol, the equivalent width ratios are 
compatible with the theoretical value of 1 .65 deduced from the 
revised atomic data, when taking into account noise features and 
the LSI 274 blend, and it is possible to obtain a reasonable si- 
multaneous fit to the two transitions. Although in the case of 
yuCol the spectrum is of excellent quality, and the two transitions 
clearly show the same structure, we could not fit them simulta- 
neously in a satisfying way (see Figure 3). More specifically, the 
narrow line in the red is correctly fitted, but the broad feature in 
the blue is not, and the probability of this being due to random 
noise is less than 5% (2cr equivalent width deficit). We suspect a 
non-interstellar origin for this broad absorption, and will discuss 
the peculiarity of this star later. 

We have also verified that for all our targets neutral titanium 
(Til) is neghgible with respect to Till (see section 1). For that 
purpose we checked for the presence of AA 3342,3636 A Til ab- 
sorption lines within the velocity interval deduced from the ob- 
served Ti II lines. We could not detect any line absorption. Since 
oscillator strengths for those Til transitions are similar to the Ti II 
ones, this implies that the contribution from Til is negligible. 

3. Comparisons between Till, DI,HI,OI, Fell 

3.1. Methods and selected targets 

In the following sections, we have combined our titanium data 
(9 targets covering logN(HI) from 19.0 to 20.98) with those of 
Prochaska et al. (7 stars with logN(HI) from 19.85 to 21.05). 



For aU abundance correlations we have adopted the D/H ratio, 
logN(HI), N(DI) and iron depletion, D(Fe), compiled in Table 1 
of Linsky et al. (2006). In addition, we have used the DI/OI ratios 
compiled in Table 12 of Oliveira et al. (2006). All the data and 
errors used in this article are listed in Table 2. For all correlations 
we have made use of the maximum amount of pubhshed data. 
Because OI and Fell measurements are not available for exactly 
the same targets, the correlation fits use different target lists. 

After completion of this work, Elhson et al. (2007) made 
available results from new Keck and UVES Till measurements. 
Some of their targets are already included in our new list. For 
those targets we have used our spectral results, resulting columns 
and errors. For those targets in common we also checked for 
consistency of the two measurements and found that our colunon 
density intervals overlap (see the caption of Table 2). There are 
4 additional stars in Ellison et al. that were not included in our 
study. We have added these 4 targets in Table 2, and we have 
performed the same correlation fits with and without those four 
targets based on the correlation analysis method discussed be- 
low. Four of our measurements are not part of the Ellison et al. 
study, the low column target aVir, the two faint subdwarfs, and 
juCol, for which Ellison et al. had an upper limit only. 

We have also compared our column density results with pre- 
vious measurements and analyses. We have disagreements with 
the results of Welsh et al. (1997). 

While for //Col we differ by 10% only, in the case of ^Pup 
our column log(TiII)= 11.18± 0.06 is a factor of 2 below their 
quoted value log(TiII)= 11.52+ 0.06. On the other hand their 
upper limits of log(TiII)= 10.79 and 10.97 for y^Vel and 0Car 
respectively are a factor of two below our measured values of 
log(Tin)=11.13±0.035 and 10.97±0.023. We note that their Nal 
and Call measurements are in excellent agreement with other 
authors (e.g. Welty, Morton and Hobbs 1996), and the discrep- 
ancy in the Till results could be due to an incomplete removal 
of the instrumental background signal which is particularly dif- 
ficult to measure at these very short wavelengths. The D/H ra- 
tio depends on the amount of molecular hydrogen. For most 
of the targets here the relative amount of molecular hydrogen 
is small and we assume N(H)=N(H1) as well as D/H= Dl/Hl. 
For three targets only (6 Car, HD90087 and PG0038+199 ) that 
have significant molecular contribution we use the numbers from 
Table 2 of Linsky et al. (2006) and use N(H)=N(HI)+2N(H2) 
and D/H= (N(DI)-hN(HD))/(N(HI)-i-2N(H2)) instead of HI and 
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Table 1. Summary of observations and titanium results. (1) distance recently revised by Millour et al (2007). *: S/N after fringe 
correction 



Target 


III 


Ml 


d(pc) 


type 


Obs-Date 


Exp-time(s) 


SN 


LogN(Till) 




J lU. i 


JU.O 




BlIII-IV 


Z.^ L/J UU 


JUA. J J 


^ou 


1 n 4.Q + n 1 ^ 


^Cen 


311.8 


1.2 


161±15 


BlIII 


15-04-06 


50x0.2 


340 


10.99 ± 0.04 


a Cm 


300.1 


-0.36 


98±6 


Bl 


10-04-06 


50x0.15 


290 


11.12 ±0.035 


7^ Vel 


262.8 


-1.1 


368(+38-13il 


WR+O 


09-04-06 


50x0.6 


320 


11.11 +0.035 


yuCol 


237.2 


-21. \ 


396±100 


09.5 


01-04-06 


40x8 


310 


11.50 ±0.02 


4" Pup 


256.0 


-4.7 


429±94 


05I 


29-03-06 


45x0.7 


250* 


11.16 ±0.06 


eCar 


289.6 


-4.9 


135±9 


BOV 


10-04-06 


40x0.9 


320* 


11.31 ±0.023 


JL9 


322.6 


-27.0 


590±160 


sdO 


19-05-06 


2x2700 


280 


11.88 ±0.10 


LS 1274 


277.0 


-5.3 


580±100 


sdB 


10-04-06 


2x2100 


110 


11.68 ±0.04 



Table 2. Compilation of D/H, D/O, titanium and iron data. Errors bars are quoted in the usual way with negative errors as indices 
and positive errors as exponants. References of Till data used in the correlations are: (a) Prochaska et al.,2005, (b), Ellison et al., 
2007, (c), this work. Targets in common with Ellison et al (2007) are marked with an asterix. The column-densities derived by these 
authors are log(TiII)= 11.10+0.02; 11.31+0.01; 11.12+0.02; 11.05+0.04; 11.07+0.04 for Vel; 9 Car; a Cru; (3 Cen and ^ Pup 
respectively and are in agreement with our determinations within uncertainties. D, H, D/H, D(Fe) are the values compiled by Linsky 
et al., (2006), while D, O, D/O are those compiled by Oliveira et al., (2006). The three targets with a very high D/O and discussed in 
sect. 3.2 are marked with a #; the first D/O value is the published value used in Figure 9a. The second number is the reduced value 
used in Fig. 9a (see section 3.2). 



Target 



log(NHI) 

17.6 

I 7 07.; 0.035 

17.884 »f 
17.93 «| 
18.06 J 

18.132g| 

18- 18 r 

18.261 

18 239 ""^^ 
18.327 «g 

18.7- 

18.76 II 

18.85 

19"' 
01 

19 11 ""25 

19.23 

19.63 SJ ! 
in 71 o.oie 
' ^ 0.0^6 

19.85 I 
19.85 Hf 
19.86 

19- 89 

19 96 " "^e 

20.03 l\ 

20.07 

20.08 
20.14»f5 

20 16 

20.16 !J f 

20 19 " "^ 

20.28 " ! 

20 4 ° 
^"•^ 0.1 

20.48 »;« 
20.78 «;«^ 
20.95 

20 98 ° 
21.05 
21.08 

21 14""^ 

'^'^■'^^ 0.06 

21 22 " 

^^■^^ 0.05 



log(DI) D/H(ppm) DI/OIxlO- D(Fe) log(Till) Till/Hl xlO'" reT 

, _ „ . — /I ni\ . , _ 1 1 _ „ n u : — rr: — n 1 1 



12.81 ii-ii^ 16:27? 



-^■^ 0.8 



450.5 

0.5 



3.5 



2.6 \i 



o q 0.6 
■^■^ 0.6 
4 0.6 
^ 0.6 
3 5 
4.2^ 
4.2 " "^ 



4.7 



3 2 
-'■^ 1.6 



7.59;3.75 
6.31;3.15 

5.62;2.81 

2 6 0.5 

358^^ 

2 5 " "^ 
"^■-^ 0.4 

2 5 
^■-^ 0.4 



-111 "1^ 

-oesSi^ 

U.UJ 

-1.07 «;| 

-0.78 «;| 
21 " 04 

-1 22 
-0.58 

I 7Q 0.05 
-'^ ■^^ 0.05 

1 7 0.04 
"'-■^ 0.04 
i 0.2 



SiriusB 
360ph 
e Eri 
31 Com 
HZ43 
Procyon 
/3Cas 
G191-B2B 
/3Gem 
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Fig. 3. 3394-3242A simultaneous line-fitting used as a consistency checked: 0Car,JL9, and LS1274. The blue-shifted absorption at 
3384A in LS1274 has no counterpart at 3242A and is very likely steUar (see text). 



DI/HI. However we have kept for figure annotations HI and 
DI/HI. 

In Figure 5 we show the Till/HI ratio as a function of the D/H 
ratio, and error bars associated with the two quantities. D/H er- 
ror bars are taken from Linsky et al. (2006). It is clear that mea- 
surement uncertainties are of the same order of magnitude for 
both related quantities. This will be the case for most attempted 
correlations. Fitting such data using standard or ordinary least 
squares methods can lead to bias in the result. In order to take 
into account both x and y errors, we have adopted the Orthogonal 
Distance Regression (ODR) method (also called "total least 
squares method") for our correlative studies. ODR minimizes the 
weighed orthogonal distance from the data to the fitted curve. In 
our case of a simple hnear fit of N data pairs ixi,yi), i=l,N with 
errors crx,- and cry,- to the function y=ax-i-b, ODR minimizes the 
quadratic sum 2((>'; - a(x; + dx/) - b)lcrylf- + (dxi/crxi)^, i.e. 
it adjusts the coefficients a and b as well as the departures dxj 
from the initial x,- positions. If the errors on x and y are identical, 
the summed quantities are proportional to the actual orthogonal 
distances between the data points and the fitted curve, hence the 
name of the method. If one of the errors is preponderant, the 
formula is equivalent to the classical least-square method. We 
have used the ODRPACK95 package (Boggs et al., 1989) imple- 
mented in the IGOR 6.0 software to perform our linear fits. In the 
case of linear fits, there is a unique solution to the minimization. 
The resulting weighted minimum distances can be visualized in 
Figure 5,8,9 and 10. 

The ODRPACK95 package works for symmetric errors, i.e. 
equal positive and negative error intervals. If there are substan- 
tial difi'erences between the positive and negative errors, the data 
point weight depends not only on the distance to the fitted line, 
but also on which side it is located with respect to it, with lower 
(resp. smaller) weight if it is located on the large (resp. small) 



uncertainty side. In order to take into account non symmetric er- 
rors we devised and we use the following iterative method: we 
first calculate the average (mean of positive and negative) error 
for Xs (resp. Ys) and run the ODR linear fit for those mean er- 
rors. A simple routine then identifies which data points are above 
or below (resp. to the left or to the right of) the fitted straight 
fine and derives from their locations which error (the positive or 
the negative one) is appropriate. A second fit is then activated 
with those new appropriate errors. This process is repeated un- 
til convergence. For linear fits and a relatively small numbers of 
points, the convergence is immediate. The method is illustrated 
in Figure 5. We show two data points for which the appropriate 
errors are of opposite signs. For all our regressions we show the 
resulting orthogonal distances between data points and the cor- 
responding adjusted locations along the fitted line. All quoted 
uncertainty intervals on the parameters are Icr and calculated 
using the error bars from the last fit after convergence. 

We list in Table 3 two diagnostics for the correlation, the 
Spearman index and the slope interval for the fitted linear rela- 
tionship. Note that the Spearman's rank correlation statistic does 
not take into account error bars, while the error on the fitted pa- 
rameters does, which means that in some cases they may give 
different results. 



3.2. Correlations of Till and Fell abundances witti D/H and 
D/0 

Visual inspection of Figure 5 shows that there is a positive corre- 
lation between Till/HI and DI/HI, with the noticeable exception 
of the star //Col. The pecuharity of this target is already visible 
in Figures 4, 6, 7 from Linsky et al. (2006) for Fell and Sill. 
For all three species the abundance is significantly above that of 
the main trend. To a lesser extent a second target also seems out 
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Fig. 4. Same as Figure 3: y^Vel,^ Pup, and ju Col. The broad feature in jx Col cannot be fitted satisfyingly (see text) 



of the main relationship, namely LS1274. In order to check the 
influence of the HI column on this result, we have displayed the 
Till columns as a function of DI in Figure 6. yuCol is clearly an 
outiier and LSI 274 is also somewhat out of the main trend. We 
have investigated possible reasons for the peculiarity of /iCol, 
but could not find anything in the literature. On the other hand, 
we have seen in the previous section that there is strong evidence 
against an interstellar origin of the broad line detected towards 
this star, on the basis of the disagreement between the 3384 and 
3242 lines. Interestingly, using the narrow line only moves the 
star into the main trend. Therefore we performed linear fits to the 
data, with and without /iCol. On the other hand we have chosen 
to maintain LS1274. Its absence does not make any noticeable 
overall difference. 

The Spearman's rank correlation statistic is 0.84 (resp. 0.59) 
without (resp. with) //Col (15 (resp. 16) points) allowing us to re- 
ject the null correlation hypothesis at better than 99.97 (resp. 98) 
%. The parameters for the slope of the linear fit are (1 .35 +0.22) 
and (1.66 ± 0.29) 10"'" respectively (see Table 3). Thus, while 
clearly /iCol modifies the fitted straight line, in both cases there 
is a positive correlation since the most probable gradient value is 
about 6 times the Icr uncertainty on this slope. The reduced 
is 1 .3 and 1 . 1 in both cases, thus there is possibly some intrinsic 
variability in addition to the positive correlation. The intercept is 
negative but compatible with zero without //Col. 

The results for Till, keeping the same methods and adding 
the 4 new Ellison et al. targets, are not noticeably diff'erent from 
those based on our own data and the Prochaska et al. data only 
(see Table 3). However, the presently derived slope is signif- 
icantly diflsrent from the one found by Ellison et al. (2007). 
There are two reasons for this disagreement, (i) We beheve their 
quoted parameters for the linear fit correspond to equal weights 
for all data points because we find the same parameters using 
their hsted data and no specific weights, (ii) The fitting method 



is not the same. Using error bars in X and Y leads to difierent 
results, (iii) The data sets are not the same. 

In order to compare in the most rehable way the Till results 
with those of Fell, we have also applied the same fitting proce- 
dure to Fell data, i.e. using the linear iron abundances (and not 
the depletion index) and the ODR method. The data and results 
are shown in Figure 7. It can be seen that the dispersion of iron 
abundances appears significantly greater than that of titanium. 
Three targets are especially clearly discrepant, with G191-B2B 
being in the most disagreement. Despite this dispersion we have 
fitted the iron abundance vs. D/H distribution without excluding 
any of the targets. We will come back to this important point 
later. Fit parameters are listed in Table 3. 

Such a positive correlation could be due to random or sys- 
tematic errors in the determination of HI (e.g. Hebrard and 
Moos, 2003) since HI appears in both related quantities and is 
by far the most difficult measurement. Indeed, changing HI may 
move the respective data point along the straight line fit. This is 
also the case for correlations between metaI/0 ratios with D/O 
ratios. Consequently, among our tests we have searched for a 
correlation between the abundance of titanium (Till/H) and the 
neutral deuterium to neutral oxygen (DI/OI) ratio. This study is 
restricted to titanium targets for which the 01 column is avail- 
able in the literature. Interestingly, although the number of data 
points used is smaller (10 points), there is clearly a trend, as 
can be seen in Figure 8. Quantitatively, the Spearman test rejects 
the absence of positive correlation at better than 98%, and the 
Till/HI vs. DI/OI slope is found to be 6.0+1.410"' (see Table 3). 
Such a trend using four independent quantities is a definite sign 
of a positive correlation between metals and deuterium. Adding 
the two Ellison et al. targets for which D/O is available does not 
significantly change this conclusion. The correlation parameters 
for the 12 point correlation are listed in Table 3. 
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Fig, 5. Ionized titanium to neutral hydrogen column ratio vs. DI/HI with (left) and without (right) the star juCol. The method used 
to take into account non symmetric error bars for the abscissa is illustrated in the left figure. For the data point 1 (resp. 2), because 
it is located on the left (resp. right) side of the fitted Une, the selected error bar is PI (resp. N2). The same method is apphed for the 
errors in the ordinate parameter. The corresponding linear fit for Fell is shown as a dotted line in the right figure. 



This result led us to correlate in a similar way the iron (Fe) 
depletion with the DI/OI ratio. Instead of using a logarithmic 
scale, we have used the linear ratio between the abundance of 
Fell and the solar abundance from Asplund et al. (2005). Here 
the number of available measurements is higher (24 lines-of- 
sight, see Table 2). 

Among these three data points there is an extreme 'outlier', 
G191-B2B. The gas towards this target is almost fully ionized, 
and as discussed in section 3.3.1 there may be a strong bias in the 
evaluation of iron depletion because iron is less rapidly ionized 
than H (see Linsky et al., 2006). We have chosen to exclude this 
target from our correlation fit, and we will come back to the ion- 
ization problem later and justify this choice. This leaves us with 
23 sight-lines, and the results are shown in Figure 9. Note that 
G191-B2B would be completely out of the figure, with a Fe/FeO 
of 0.26. As seen with titanium in Figure 7, the data point distri- 
bution reveals a clear positive correlation. Three points show a 
rather large deviation from the mean curve and seem to have a 
much higher D/0 compared to the other targets. Those very dis- 
tant targets are indeed the only three targets for which very high 
D/O ratios have been claimed: WD 1034-1-001, BD -1-39 3226 and 
TDl 32709 (Oliveira et al. 2006). However, toward those targets 
the /1974A OI transition has not been used in the OI colunon 
determination, or was not significantly detected. Not using this 
transition could lead to underestimation of N(OI) by factors of 
^2, as it has been shown in the cases of Feige 110 (Hebrard et al. 
2005) and LSE 44 (Friedman et al. 2006). Using this transition, 
the OI column density toward BD -h39 3226 is also multiplied 
by a factor of 2.6 (Oliveira et al. 2006). Thus, as suggested by 
Figure 9, the D/O measurements might be overestimated by a 
factor of 2 to 3 for these three hues of sights. 



We have thus performed two analyses: firstly using the 3 high 
Dl/01 values from the strong line fitting, and secondly after divi- 
sion by a factor of 2 of the DI/OI ratios for only these three dis- 
tant stars, hypothesizing that a correction similar to the Feige 1 10 
re-evaluation is likely to be necessary. When this allowance for 
correction of line saturation eff'ects is included, the three corre- 
sponding data points all move back closer to the main stream of 
data points. We have also performed a similar linear fit, exclud- 
ing the three targets instead of revising their DI/OI and found 
results very similar to those we show here. The two fits are dis- 
played in Figure 9 and the linear fit parameters before and af- 
ter correction are given in Table 3. The Spearman test rejects 
the absence of correlation at better than 99.6% and 99.99 % in 
the first and second case respectively. The gradient is somewhat 
steeper after the correction for the three high colunms, namely 
1.49 ±0.26 instead of 1.28 ±0.22. 



The Till and Fell positive correlations with D/O demonstrate 
three important points. Firstly, potential errors in the measured 
values of HI as the source of a link between D and metals in the 
gas phase is a hypothesis that can now be excluded as the main 
origin of the deuterium-metal correlation. If errors in the mea- 
surement of HI were the case, we would simply have randomly 
distributed points. Secondly, the D/O ratio, despite its difi^erent 
evolution as a function of the gas column compared to DI/HI, is 
shown to vary in a specific way, and shows a link to the metal 
abundance. Third, D/O cannot be linked to metal depletion prin- 
cipally through astration effects, since this would produce the 
opposite trend, with D/O being anti-correlated with the degree 
of astration and thus with metal abundance. 
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Table 3. Correlation parameters. LHW= this paper, PTH05 = Prochaska et al. (2005), EPL07= Ellison et al. (2007), Lin06 = Linsky 

et al. (2006), O106 = Oliveira et al., (2006) 
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Fig. 7. Ionized iron abundance normalised to (Fe/H)o=-4.55 
(Asplund et al., 2005) vs. DI/HI. This figure is identical to Figure 

3 from Linsky et al.(2006) except for the linear scale. The best 
fit linear correlation using the Orthogonal Distance Regression 
(ODR) is shown as a solid line, and for comparison the titanium 
gradient obtained with the same method (Figure 5) is shown as 
a dashed line after the appropriate scaling to obtain a similar 
ordinate for D/H=20. It can be seen that the dispersion is sig- 
nificantly larger for Fell than for Till. Lines-of-sight filled with 
strongly ionized gas (G191-B2B, 36 Oph) are positive extreme 
"outliers", while those corresponding to very large colunms tend 
to produce negative "outliers". 



3.3. Gradient comparisons 

3.3.1 . Comparisons between titanium and iron 

It is interesting to compare the variation amplitude of titanium 
and iron. Titanium is strongly refractory and one expects a 
higher amplitude of the gas-phase abundance variation, as ar- 



gued by Ellison et al. (2007). Figure 5 shows the Fell vs. DI/HI 
gradient from Figure 7 superimposed on the Till vs. DI/HI gra- 
dient. Contrary to expectations based on condensation temper- 
atures. Till varies slightly less than Fell, although as shown in 
Table 3 and taking into account the two error bars, the two gra- 
dient intervals nearly overlap. In the case of Dl/Ol the two gra- 
dients are compatible, but the Till vs. DI/OI gradient is less pre- 
cisely defined than for DI/HI. 

We believe that the absence of a stronger gradient for Till is 
not an inevitable concern with respect to the deuterium depletion 
hypothesis of Linsky et al (2006) and that the observed superior- 
ity or similarity of the iron slope can be explained by ionization 
effects. Linsky et al. (2006) have already discussed the potential 
impact of diflFerential ionization of Fell and HI on the Fell/Dl ra- 
tio variation. The ionization potential of Fell is 16.2 eV, and as a 
consequence hydrogen may ecome ionized by the radiation field 
while Fell remains only slightly affected. The case of G191-B2B 
, the extreme outlier, is typical. From their analysis of high res- 
olution HST-GHRS spectra, Lemoine et al (1999) inferred the 
existence of two main clouds along the LOS to the hot white 
dwarf. For the main HI cloud (B) the Fell/HI ratio is 0.16 lO'^ 
while for the second, strongly ionized cloud (A), the same ratio 
is 2.25 IQ-^, i.e. 13-14 times above cloud B. The existence of the 
ionized fraction contained in cloud A has the effect of increasing 
the overall ratio (for A+B) by a factor of more than 3 compared 
to cloud A alone. This is a very strong effect, and it is clear from 
Figure 7 that Unes-of-sight filled with significantly ionized gas, 
such as G191-B2B or 36 Oph have a systematically higher singly 
ionized iron abundance, that probably does not reflect the actual 
iron abundance. As a consequence it is likely that a high ion- 
ization degree, a phenomenon expected to be present in shock- 
and photo-heated regions, i.e. in regions with the strongest evap- 
oration and deuterium release from grains, also tends to increase 
the Fell/HI ratio. In other words, in the absence of precise cor- 
rections for the gas ionized fraction, iron release in the gaseous 
phase and apparent abundance increase due to ionization bal- 
ance effects act in a cumulative way. This tends to increase the 
Fell/H ratio and the larger gradient for Fell may result from this 
trend. As we explained earUer G191-B2B is an extreme case and 
we have removed it from the correlation with DI/OI shown in 
Figure 9, which has the effect of decreasing the slope. Probably 
other remaining targets such as WD221 1-495 have an effect on 
the slope of the linear fit. G191B2B is also the most discrepant 
point for the Fell/HI vs. DI/HI correlation, but because there are 
two times more targets for this correlation it is not isolated, and 
this is why we have kept it. Looking at Figure 7 it is clear that 
three data points including this star tend to increase the gradi- 
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Fig. 8. Left: Ionized titanium vs. DI/OI. The corresponding lin- 
ear fit for Fell (Figure 5) is shown for comparison. 



ent of the correlation. These ionization effects are absent in the 
case of Till, due to the similarity of the potentials. This is why 
we believe that the absence of the expected stronger gradient for 
titanium compared to iron is due to an overestimate of the latter. 

3.3.2. Comparison between the titanium and iron gradients 

vs. D/H and vs. D/0 

While the positive correlation between the metal abundance and 
D/O can not be explained by astration (it would be a negative 
trend), it may be well explained, as for D/H (Linsky et al., 2006), 
by the depletion hypothesis. Here we compare the two gradients 
obtained from the previous correlations of Till and Fell with 
D/H and D/O. 

Along with deuterium and metals that evaporate from grains 
and are released in the gaseous phase, it is reasonable to expect 
some simultaneous release of oxygen. However the relative in- 
crease of the 01 concentration due to this release is much smaller 
than for Dl, because the oxygen abundance is the gaseous phase 
is much higher than for deuterium. Also, data show a quasi- 
constancy of the OI/HI ratio, showing that an increase due to 
dust evaporation is negligible. Therefore, in the depletion sce- 
nario, for a given increase of the metal gaseous abundance, i.e., 
of dust evaporation, DI/OI should increase in the same propor- 
tions as DI/HI. Figures 5,7,8, and 9 and Table 3 show that this is 
the case globally, i.e. the metal abundance, the D/H ratio and the 
D/O ratio all vary by a factor of 3-3.5, and there are significant 
linear correlations between metals and D/H on one hand, metals 
and D/O on the other. However, there is no simple proportional- 
ity between the three quantities, since some lines of sight with 
very high D/H ratios do not present correspondingly high D/O 
values (Hebrard et al, 2005). 



It is striking that error intervals for the fit parameters allow 
(or nearly allow) a simple proportionality between Ti or Fe and 
D ( the linear fit in this case goes through zero). This implies 
that D increases as rapidly as metals. In the framework of the 
deuterium depletion scenario of Draine (2004) and Linsky et 
al. (2006), one would expect a preferential release of deuterium 
from PAHs and the carbonaceous envelopes of grains, whereas 
metals (O and Fe) are instead preferentially removed from the 
cores of interstellar grains, and therefore a strict proportionality 
is somewhat surprising. 

The comparison between the residuals for Till/HI vs. D/H 
(X^/N = 1.3 and 1.1, see Table 3) and D/O (x^/N = 2.7 and 2.3) 
shows that the fit is of significantly better quality vs. D/H than 
vs. D/O. We believe this is due to the fact that in the former 
case N(H) appears in the two correlated quantities. As a con- 
sequence, measurement errors on H tend to produce dispersion 
along the fitted line, that does not increase the x^- This shows 
that a fraction of the correlation between Till/Hl and D/H is due 
to uncertainties on HI columns. In the latter case of the correla- 
tion with D/O, measurement errors in H bring a dispersion along 
the ordinate axis that appear in the x^- Such effects for iron are 
not seen because the residuals are dominated by the ionisation 
biases that introduce a strong dispersion. Especially, we have re- 
moved G191B2B for the correlation with D/O and kept it for 
D/H, which explains the higher x^/N for D/H. 

3.3.3. Correlation between titanium abundances and 
column-densities 

We have searched for deviations from proportionality and for a 
potential link between the titanium to deuterium abundance ratio 
and the HI column. Linsky et al.(2006) have showed evidence for 
an increase of the DI/Fell ratio along with the HI column. More 
precisely, there seems to be an increase from very low columns 
(10'^ cm"^) to columns of the order of lo'^^'^-^ cm"^ and con- 
stancy above log(NH) ^ 19.5 cm"^. They interpreted this trend 
as evidence for a quasi proportionality of D and Fe abundances 
in general, and attributed the low to medium column increase as 
being due to significant ionisation of the diffuse gas close to the 
Sun, an effect similar to what we have invoked for the Ti vs. Fe 
comparison. As discussed by the authors and previously in this 
paper for HI, in ionized gas D may be more easily ionised than 
Fell, leading to an overestimate of the abundance of Fell at low 
columns within the Local Bubble. 

In the case of Till, ionisation effects are not supposed to have 
any contribution due to the coincidence between the two ionisa- 
tion potentials of Till and HI (or DI). It is therefore particularly 
interesting to perform the same correlation test as the one done 
by Linsky et al. but for Till instead of Fell. Figure 10 shows the 
Till/DI ratio as a function of logN(HI). Similarly to the Linsky 
et al. (2006) results for Fell, the smallest HI columns seem to 
correspond to the highest Tin /DI ratio. While the trend is not 
as clearly defined as in the correlation analyses in the previous 
sections, with a Spearman rejection of null correlation at about 
75% probability only (and a negative result if we do not include 
the high column targets of Elhson et al (2007), see Table 3), the 
least-squares Unear fit with the ODR method taking into account 
error intervals on TII/DI and logN(HI) demonstrates a negative 
slope at better than 3.5 crs. However, the slope is very strongly 
dependent on the three points with high columns and small error 
bars and more or better data are needed to assess this result. 

The sign of this inferred Till/DI gradient cannot be explained 
by less astration for high columns, i.e. far from the Sun, since it 
would contradict the increase of OI/HI at large columns found by 
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Fig. 9. Iron gaseous abundance (normalised to the Asplund et al abundance) and Deuterium to Oxygen ratio, a: Published DI/OI 
results, b: After division by a factor of 2 of the three DI/OI ratios for the three distant stars (see text). 
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Fig. 10. Till/DI ratio as a function of the HI column: all available measurements. The linear fit using the Orthogonal Distance 
Regression method is shown as a solid line. The slopes found by Linsky et al. (2006) for the Fell/DI ratio are shown as dashed and 
dotted lines (resp. all data, and data with logN(HI) > lO'^ cm"^). 



Andre et al. (2003) and more recently by Oliveira et al. (2006). 
Tentatively, the sign of the slope may correspond to preferential 
evaporation of deuterium compared to titanium at the highest 
columns, i.e. in very dense and colder regions which become 
preponderant along very distant lines-of-sight. In other words, D 
in this case is released into the ISM before metals, and only in 
warmer and more tenuous media are the release rates identical. 

Superimposed on the Tii/DI vs. logN(III) fitted relationship 
are shown the corresponding results of Linsky et al. (2006) for 
Fell. Ordinates have been multiplied by a constant, and we are 
interested here in the slope only. Our results that correspond to 



logN(HI) > 19 should be compared to the flat curve calculated 
by these authors for targets with logN(HI) > 19.3. 



As we already mentioned, titanium data are very appropriate 
here because of the absence of ionization effects that may in- 
troduce biases. However there is a lack of measurements with 
logN(H) smaller than 10^^ cm"^ to confirm unambiguously a 
trend. 
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4. Discussion and Conclusions 

We have observed nine southern hemisphere D/H target stars at 
high spectral resolution and derived the column-densities of in- 
terstellar ionised titanium towards those targets. A comparison 
between the D/H ratio and the titanium depletion reveals a clear 
correlation and reinforces the result of Prochaska et al. (2005) 
and the more recent result of Ellison et al. (2007). Our work 
extends the correlation over a broader range of interstellar HI 
columns, namely more than two orders of magnitude, showing 
that the trend extends down to log N(HI)=19.0 cm"^. We also 
make use of the orthogonal distance regression method for our 
hnear fits, in order to take into account errors in quantities both 
in abscissa and ordinate. 

One star is far out of the correlation, yuCol. This target is 
already strongly discrepant in correlations shown by Linsky et al 
(2006) between D/H and iron and silicon abundances. We have 
shown evidence here for a non-intersteUar origin of the broadest 
absorption feature and note that excluding this absorption would 
bring the star back lo the main trend. 

We have compared both the Ti and Fe depletions with the 
DI/HI ratio, but also with the DI/OI ratio, when available. For 
both metals there is a clear relationship also with DI/OI, which 
demonstrates that the metal abundance- D/H correlation is not 
primarily due to uncertainties in HI column measurements, and 
that D/O is also correlated with metal depletions in the ISM. 

The ionized titanium correlation gradient is found to be simi- 
lar to or slightly smaller than the corresponding gradient for ion- 
ized iron. Ellison et al (2007) have emphasized that given its high 
condensation temperature one would expect the titanium gradi- 
ent to be the highest. We have argued here that ionization effects 
do affect Fell/HI and result in an overestimate of the gradient, 
while the titanium abundance is unaffected due to the nearly ex- 
act coincidence between Tin ad HI ionization potentials. Thus 
we beUeve that the slopes measured for both metals are not a 
concern with respect to the depletion hypothesis as the source 
of deuterium abundance variations. Further work is needed to 
correct for the ionization and obtain a more reliable comparison 
between the different metal-deuterium relationships. 

There is some (although weak) evidence that the Till/DI ra- 
tio decreases towards high column densities. This may signify 
that D is preferentially released into the ISM (when compared to 
titanium) in the dense regions that are predominant along distant 
Unes-of-sight. It is important to investigate further this potential 
relationship, because, if confirmed, it would be a first indication 
of the different sources of interstellar D and metals. Recording 
more titanium data along shorter lines-of-sight is essential. TUI, 
that does not suffer ionisation biases, is an ideal tracer here. 

It remains to explain why very high D/H ratios do not have 
corresponding very high D/Os, as already argued by Hebrard et 
al (2005). This is reinforced here by our findings (section 3.2 
and Figure 9a) that the three high D/O values found by Oliveira 
et al (2006) may be overestimated. Additionally, Oliveira and 
Hebrard (2006) derive a high D/H for a very distant sighthne 
while one would expect in this case a very strong depletion. 
Certainly some of the physical processes are not yet understood. 
On the other hand, since very discrepant points appear in the cor- 
relation with D/H, but none with D/O, the high D/H ratios may 
be a consequence of underestimated H columns. 

Assuming proportionality between metals and D depletion, 
one can derive a mean O/H ratio from the ratio of the Ti/H vs. 
D/H and Ti/H vs. D/O slopes (or similarly for Fe). Using Till one 
finds (Table 3) 6.05 10"Vl-35 lO""* = 448 ppm using our targets 
and the Prochaska et al. targets, and 6.46 10 Vl-29 10 = 501 



ppm using all targets. Using Fell one finds 1.49/4.96 10^= 300 
ppm. There must be some significance to this difference that re- 
mains to be elucidated. These two values are respectively higher 
and lower than the mean value of 347 ppm derived by Cartledge 
et al. (2004) within 800 parsecs. It is also interesting to compare 
with the range of solar photospheric O/H values recently pro- 
posed by Grevesse et al. (2007). The value derived here from the 
Till data, 480-500 ppm, is within their quoted range of 407-513 
ppm, while the value derived from Fell, 300 ppm, is significantly 
below. These calculations suggest that the slopes derived from 
the Uneiir fits using Till are meaningful. 

Is there something to learn about the local ISM history from 
the D/H abundance pattern? The large-scale event that has given 
rise to the Gould belt of enhanced star formation is about 60 
Myrs old (see e.g. Perrot & Grenier, 2003). Two main classes 
of scenarios have been proposed for the belt origin: the inter- 
action between a giant external cloud and the disk on one hand 
(Comeron and Tora, 1994, Olano, 2001) and a strong explosive 
event on the other (Olano, 1982). Because the age of the belt is 
smaller than the mixing time of gases (a few hundreds Myrs, see 
e.g. simulations by De Avillez and Breitschwerdt, 2004), in the 
former case of cloud collision the belt may have left imprints in 
the form of incomplete mixing between disk and "fresh" gas. In 
the latter case there must be dust evaporation at the shocked front 
regions of the belt but no gas mixing. A link between the Gould 
belt expansion and the D/H distribution has been suggested be- 
cause highly variable D/H regions shown by Linsky et al (2006, 
Fig. 1) are located at distances (inferred from the H columns) 
that roughly correspond to the front of the expanding Gould belt 
of young stars and supernovae (Lallement, 2007). Whatever the 
origin of the belt, the subsequent and observed expansion pro- 
duces star formation and supernovae at the periphery of the belt, 
i.e. evaporation of dust in the shocked and heated regions of its 
periphery, which could explain the observed variabiUty of D/H 
at intermediated columns (^iQ'^^-^-^^.s ^jj^-z-j ^he difference be- 
tween the two scenarios for the belt origin is that in the cloud 
collision case some deuterium abundance variabihty due to gas 
mixing must add to abundance variation due to grain destruction, 
and there must be a larger amount of less astrated gas within the 
Gould belt as compared with outside the belt. There are no signs 
of astration variations that follow from the present correlative 
analysis, in agreement with the constancy of O/H up to large 
distances (Andre et al., 2003). If D and metal abundance vari- 
ability has some link with the belt, then the absence of signs of 
differential astration at columns lower than 10^' cm"'^ seems to 
favor the explosive event scenario. 

Independantly of a potential effect of the Gould belt on the 
D/H distribution, our results provide a support for the ideas, orig- 
inally suggested by Draine et al. (2004) and Linsky et al. (2006), 
that time-dependent deuterium depletion onto dust grains is re- 
sponsible for the variabihty of the D/H ratio, and for the low 
values of D/H observed for many lines of sight beyond the Local 
Bubble. On the other hand, it remains to explain why very high 
D/H ratios do not correspondingly present very high D/O values. 
This is important for the precise determination of the present 
galactic D/H ratio. 
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